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Abstract — This paper considers an orthogonal frequency divi- 
sion multiplexing (OFDM) point-to-point wireless communication 
system with simultaneous wireless information and power transfer. 
We study a receiver which is able to harvest energy from 
the desired signal, noise, and interference. In particular, we 
consider a power splitting receiver which dynamically splits the 
received power into two power streams for information decoding 
and energy harvesting. We design power allocation algorithms 
maximizing the spectral efficiency (bit/s/Hz) of data transmission. 
In particular, the algorithm design is formulated as a nonconvex 
optimization problem which takes into account the constraint 
on the minimum power delivered to the receiver. The problem 
is solved by using convex optimization techniques and a one- 
dimensional search. The optimal power allocation algorithm serves 
as a system benchmark scheme due to its high complexity. To strike 
a balance between system performance and computational com- 
plexity, we also propose two suboptimal algorithms which require 
a low computational complexity. Simulation results demonstrate 
the excellent performance of the proposed suboptimal algorithms. 

I. Introduction 

Orthogonal frequency division multiplexing (OFDM) is a 
promising air interface to fulfill the growing demands for a high 
spectral efficiency, due to its flexibility in resource allocation 
and high resistance against channel delay spread. Nevertheless, 
in energy limited wireless networks, the lifetime of communi- 
cation nodes remains the bottleneck in guaranteeing quality of 
service (QoS) due to the constrained energy supply. Recently, 
energy harvesting technology has received considerable interest 
from both industry and academia. In particular, it has become 
a viable solution for prolonging the lifetime of networks since 
it provides a perpetual energy source and self-sustainability to 
systems Q]-|6). 

Traditionally, tides, solar, and wind are the major natural 
renewable energy sources for energy harvesting. Unfortunately, 
the availability of these energy sources is limited by climate or 
location which may be problematic in indoor environments. On 
the other hand, harvesting energy from ambient radio signals in 
radio frequency (RF) introduces a new paradigm of energy man- 
agement. More importantly, wireless energy harvesting provides 
the possibility for simultaneous wireless information and power 
transfer |3|-|6|. Nevertheless, this combination poses many new 
challenges for the design of resource allocation algorithms 
and receivers. In J5] and J4], the optimal tradeoff between 
power and information transfer was studied for different system 
settings. However, the receivers in and J4] are assumed 
to be able to decode information and extract power from the 
same received signal, which cannot be achieved in practice 
yet. Consequently, a power splitting receiver was proposed in 
(5) and J6) for facilitating simultaneous information decoding 
and energy harvesting. In particular, the authors in studied 
the resource allocation algorithm design for power splitting 
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receivers of single carrier systems in ergodic fading channels. 
Yet, the assumption of channel ergodicity may not be justified 
in slow fading channels and the results in J6) may not be 
applicable to multicarrier systems. Besides, an algorithm for 
maximizing the spectral efficiency of a system with power 
splitting receiver has not been reported in the literature so far. 

Motivated by the aforementioned prior works, in this paper, 
we first derive an optimal algorithm for maximizing the system 
spectral efficiency in slow fading channels. Due to the associ- 
ated high computational complexity of the optimal algorithm, 
we propose two suboptimal resource allocation algorithms 
with low computational complexity which are based on the 
coordinate ascent method and convex optimization techniques. 

II. System Model 
In this section, we present the adopted system model. 

A. OFDM Channel Model 

A point-to-point OFDM system consisting of one transmitter 
and one receiver is considered. We assume that the receiver 
is able to decode information and harvest energy from noise 
and radio signals (desired signal and interference signal). All 
transceivers are equipped with a single antenna, cf. Figure Q] 
There are rip subcarriers sharing the system bandwidth of B 
Hertz; each subcarrier has a bandwidth of W = B/uf Hertz. 
The channel impulse response is assumed to be time invariant 
(slow fading) and the channel gain is available at the transmitter 
for resource allocation purpose. In addition, the receiver is 
impaired by a co-channel interference signal emitted by an 
unintended transmitter. The downlink received symbol at the 
receiver on subcarrier i G {1, . . . , uf} is given by 



Yi = VP'giHtXt +h + Z\ + Zf, 



(1) 



where Xi, Pi, and Hi are the transmitted symbol, the transmit- 
ted power, and the multipath fading coefficient on subcarrier 
i, respectively. I and g denote the path loss attenuation and 
shadowing between the transmitter and receiver, respectively. 
Zf is the signal processing noise on subcarrier i with zero mean 
and variance a\ s . Zf is caused by quantization errors. Zf is 
the antenna noise on subcarrier i which is modeled as additive 
white Gaussian noise (AWGN) with zero mean and variance 
a^a. It is the received interference signal on subcarrier i with 
zero mean and variance cr| . 

B. Hybrid Energy Harvesting and Information Receiver 

In practice, the signal used for information decoding cannot 
be reused for harvesting energy due to hardware limitations 
J6). As a result, we follow a similar approach as in ||6] and 
focus on a power splitting receiver to facilitate the concurrent 
information decoding and energy harvesting. In particular, the 
receiver splits the received signal into two power streams 
carrying proportions of pi and 1 — pi of the total received 
signal power before any active analog/digial signal processing is 
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Fig. 1. Communication system model for simultaneous wireless information and power transfer. 



performed, cf. Figure Q] Consequently, the two streams carrying 
a fraction of pi and 1 — pi of the total received signal power 
are used for decoding the information embedded in the signal 
and energy harvesting, respectively^. 

III. Resource Allocation 

In this section, we first introduce the system capacity and 
then formulate the corresponding resource allocation algorithm 
design as optimization problem. 

A. Instantaneous Channel Capacity 

The channel capacity (maximum spectral efficiency) between 
the transmitter and the receiver on subcarrier i with channel 
bandwidth W is given by 

(2) 
(3) 



C 2 = Wlog 2 [1 + PiTi) with 
p l h\H l \ 2 



r 



Pi (a 2 z a + erf.) + ol 



where P;T; is the received signal-to-interference-plus-noise 
ratio (SINR) on subcarrier i. Here, the interference and signal 
processing noise are treated as AWGN for simplifying the al- 
gorithm design. The system capacity is defined as the aggregate 
number of bits delivered to the receiver over n p subcarriers and 
is given by 



11 f 

i=l 



(4) 



where V = {Pi > 0,Vi} is the power allocation policy and 
p = {0 < pi < 1, Vi} is the power splitting ratio policy. 

B. Optimization Problem Formulation 

The optimal resource allocation policy, {V*, p*}, can be 
obtained by solving the following optimization problem: 



max U(V,p) 



up 



s.t. c\:Yp l \H l \ 2 gHi- P i)>P r ^ 



2 = 1 

n p 



C2: J2 p i < p max, C3:P C + J2 eP * ^ P pg , 



i=l 



C4: Pi = p h \jj^i, C5: P > 0, V*, 
C6: < pi < 1, Vi. 



(5) 



Here, P mi „ in CI is the minimum required power transfer to 
the receiver which represents a QoS requirement, -q denotes 

'For the sake of presentation, we assume that the power splitting ratio can 
be different across different subcarriers at the moment. The implementation 
constraints on the power splitting will be taken into account when we introduce 
the problem formulation for resource allocation algorithm design. 



the energy harvesting efficiency of the receiver in converting 
the received radio signal to electrical energy for storage. C2 
constrains the maximum transmit power of the transmitter such 
that it will not exceed P max - In practice, the value of P max 
is related to hardware limitations of the transmitter and/or the 
maximum spectrum mask specified by regulations. Variables 
Pc and e in C3 are two constants account for the circuit 
power consumption in the transmitter and the inefficiency of the 
power amplifier, respectively. C3 indicates that the maximum 
power supply from the power grid is Ppc, cf. Figure Q] and 
the total power consumption of the transmitter is controlled to 
be less than Ppc- C4 accounts for the hardware limitations 
of the power splitting receiver. In particular, pi is required to 
be identical for all subcarriers. Otherwise, an analog adaptive 
passive frequency selective power splitter is required at the 
receiver which results in a high system complexity. 

IV. Solution of the Optimization Problem 

Problem <j3j belongs to the class of nonconvex optimization 
problems. In particular, the power splitting ratio, pi, couples 
with the power allocation variable, Pi, in the SINR on each 
subcarrier which makes (fj) a nonconvex function with respect 
to (w.r.t.) pi and Pi. As a result, efficient convex optimization 
techniques may not be applicable for obtaining the global 
optimal solution. In the following, we propose an optimal 
algorithm and two suboptimal algorithms for system capacity 
maximization. The optimal resource allocation algorithm com- 
prises a full search for pi and convex optimization techniques. 
Specifically, we maximize the system capacity w.r.t. the transmit 
power for a given fixed pi. Then, we repeat the procedure for 
all possible values of pi and record the corresponding achieved 
system capacitiefl At the end, we select that pi as the optimal 
power splitting ratio from all the trials which provides the 
maximum system capacity. We note that although the optimal 
resource allocation algorithm achieves the global optimal sys- 
tem performance, it incurs a prohibitively high computational 
complexity to the transmitter which is not desirable for time 
constrained wireless communication services. 

A. Optimal Algorithm 

In this subsection, we solve the power allocation optimization 
problem by convex optimization techniques for a given set of 
p. To this end, we first obtain the Lagrangian function of (0: 

£(A,/?, 7 ,P) (6) 
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2 In general, an np dimensions full search is required to obtain the optimal 
power splitting ratio. Yet, the search space can be reduced to a one-dimensional 
search because of constraint C4. 



TABLE I 

Iterative Resource Allocation Algorithm. 



Algorithm Suboptimal Resource Allocation 1 
1: Initialization: N max = the maximum number of iterations 

and A = the maximum tolerance 
2: Set iteration index n = and initial resource allocation 

policy {V n ,p n } 
3: repeat {Iteration} 

4: For a given set of p n , obtain an intermediate power 

allocation V' n from (O 
5: For a given set of V' n , obtain an intermediate power 

splitting ratio p' n from (© 
6: if \U{V' n , P ' n )-U{V n , Pn )\< A then 
7: Convergence = true 
8: return {V*, p*} = {V, p' n ) 
9: else 

10: Set {V n , p n } = {V' n , p'J and n = n + 1 
li: Convergence = false 
12: end if 

13: until Convergence = true or n = -/V max 

Here, 7, f3, A > are the Lagrange multipliers associated 
with the transmitter power usage constraints CI, C2, and C3, 
respectively. The non-negative transmit power constraint in C5 
will be captured in the Karush-Kuhn-Tucker (KKT) condition^ 
when we derive the optimal transmit power. We note that 
constraints C4 and C6 do not contribute to the Lagrangian 
function since they are independent of Pi. Yet, these constraints 
will be considered in the full search over pi. 

Then, by using the KKT conditions for a fixed set of 
Lagrange multipliers, the optimal power allocation on subcarrier 
i is obtained as 



P* = 



W 



ln(2)(a + /3- 7 ) 



Pilg\Hi 



V*, (7) 



and [x] + = max{0, x}. It can be observed from © that the 
power allocation solution is in the form of water-filling. Since 
Pi is fixed in this optimization framework, the influence of pi is 
treated as part of the channel gain. Lagrange multiplier 7 forces 
the transmitter to allocate more power for transmission to fulfill 
the minimum power transfer requirement P[^ r The optimal 
Lagrange multipliers can be easily found via the gradient 
method or off-the-shelf numerical solvers |0, JS]. 

In the following, we focus on the design of two suboptimal 
power allocation algorithms with low computational complexity 
compared to the optimal power allocation algorithm. 

B. Suboptimal Algorithm 1 

The problem formulation in d5} is concave w.r.t. V or p, but 
not jointly concave w.r.t. both of them. As a result, an iterative 
coordinate ascent method is proposed to obtain a locally optimal 
solution of (0 and the algorithm is summarized in Table 
U V n and p n denote the power allocation and power splitting 
policy in the nth iteration, respectively. The overall algorithm is 
implemented by a loop which solves two optimization problems 
iteratively. In each iteration, we execute line 4 in Table [Qby first 
keeping the power splitting ratio p n fixed and optimizing V n via 
10. Then in line 5, we keep the updated transmit power V n fixed 
and optimize p n . The procedure iterates between line 4 and 

3 Since the problem is concave w.r.t. V for a given set of p and it satisfies 
Slater's constraint qualification, the KKT conditions provide necessary and 
sufficient conditions for the optimal transmit power allocation. 



line 5 until the algorithm converges or the maximum number 
of iterations has been reached. We note that the convergence to 
a local optimal solution is guaranteed Q. 

For solving the power splitting ratio p n with a fixed V n in 
line 5, we apply the KKT conditions for 10 which yields 



9,/2/(^+ g f i )/(7 + E^E,a j ) 
v/E(2) {(a^+aD + mHgP^ 



la 2 + lW,P,o» 



{0% 



H 



n g Pi) 



(8) 



where 9, = ^H^lgP^s (7 + £ 4 = 

AW(a 2 a +olf+AW\H l \ 2 lgP l {o- 2 za +al ) + \H l \ 2 lgP l a 2 s (7+ 
2~2jjLi 2~2i d,j) m (2), and is a two-dimensional Lagrange 
multiplier chosen to satisfy the consensus constraint C4. Opera- 
tor [x] , is defined as [x] . = c, if x > c, [x] , = x, if d < x < 
c, [x] , = d, if d > x, respectively. It can be verified that p* is a 

monotonic increasing function of Pi, i.e., ^p- > 0. As a result, 
we expect that when Pi is large enough, p* will have a value 
close to 1. We note that the proposed suboptimal algorithm 
has a polynomial time complexity due to the convexity of the 
problem formulation w.r.t. the individual optimization variables. 

C. Suboptimal Algorithm 2 

In this section, we propose a suboptimal resource allocation 
algorithm which is asymptotically optimal in the high SINR 
regime. For facilitating the design of an efficient resource 
allocation algorithm, we augment the optimization variable 
space by replacing pi with two auxiliary variables, i.e., p\ 
and pf . Specifically, p\ and pf are associated with the power 
splitting ratios of the power streams for information decoding 
and energy harvesting, respectively. Note that p\ + pf = 1 
has to be satisfied which indicates that the power splitting unit 
does not introduce any extra power gain to the received signal. 
In addition, we rewrite and approximate the channel capacity 
between the transmitter and the receiver on subcarrier i as 

PiPligm 2 



Ci 



Wlog 2 ( 



1 



PiPligm 2 



(9) 



in high SINR, i.e., log 2 (l + x) « log 2 (x),x » 1. As a 
result, the objective function is now jointly concave w.r.t. p\ 
and Pi since the two eigenvalues of the Hessian matrix of Ci, 

-1 ari d (^.) a +pi(2^+2<K. 



-3-, are non-positive. 



P, log(2) (p 7)2 log(2 )^2 3+p 7 a 2 a+p I a 2^ 

Furthermore, we rewrite constraints CI, C4, and C6 as 



preq 
min 

Pf ' 



CI :J2 p i\ H i\ 2 9 l V> 

»=i 

C4:pf =p?,\Jr = {2,...,n F }, 
C6:p\ +pf = l,Vi, 



and 



(10) 

(ID 
(12) 



respectively. Finally, we impose a non-negative value constraint 
on the auxiliary variables, C7: pf,pf > 0,Vz. Hence, the con- 
straints C1-C7 span a convex feasible solution set. Therefore, 
the problem with the approximated channel capacity ((9) is 
jointly concave w.r.t. the optimization variables and traditional 
convex optimization techniques can be used for obtaining the 



TABLE II 

System parameters 



Receiver distance 


10 meters 


Multipath fading distribution 


Rician fading (Rician factor 6 dB) 


Carrier center frequency 


470 MHz 


Number of subcarriers np 


128 


Total bandwidth 


20 MHz 


Sisnal processing noise a 2, s 
- - S 


-35 dBm 


Antenna noise crf a 


-115 dBm 


Channel path loss model 


TGn path loss model 


Lognormal shadowing g 


1 


Circuit power consumption Pq 


40 dBm 


Max. power grid supply Pq 


50 dBm 


Power amplifier power efficiency 


1/e = 0.16 


Min. required power transfer P' "f 


dBm 


Energy harvesting efficiency rj 


0.8 



solution. The Lagrangian of the transformed problem with the 
approximated objective function is given by 



= £ Q - A (Pc + X; ^ - Pre) - J2 (pf - Pr) 

i=l i=l r=2 

j=i i=i 

np 

-X>(pf+pf-i), ( 13 > 



2 = 1 



where and 4> are the Lagrange multiplier vectors with 
elements pi and <j>j which are associated with constraints C6 
and C4, respectively. Therefore, the transmit power and the 
power splitting factor can be obtained as 



W 



ln(2)(a + /3- 7 ) 



Vi, 



Mi 



{AW{&+o* u )+ko*. ln(2)) 



2^(^+4)71^(2) 



, Vi, 



(14) 



(15) 



pf* 



2(o*.+of 4 )_ 

, Vr{2, . . . ,Uf}, and 



P 



req 



Mi 



.7 = 1 ' 



(16) 



(17) 



The power allocation solution in ( fT~4-b suggests that equal power 
allocation across different subcarriers is optimal in the high 
SINR regime. Moreover, P* is decoupled from pf* and, in 
contrast to suboptimal algorithm 1, iteration between P* and 
p\ * is not required for obtaining the solution. Since the consid- 
ered problem with the approximated objective function is jointly 
concave w.r.t. the optimization variables, the problem can be 
solved efficiently by finding the optimal Lagrange multipliers 
with numerical solvers J8]. 

V. Simulations 

In this section, we evaluate the performance of the proposed 
power allocation algorithms using simulations. The simulation 
parameters can be found in Table [TT] For the optimal resource 
allocation algorithm, we use 1000 equally spaced intervals 
for quantizing the range of pi for facilitating the full search. 
Besides, the maximum number of iterations iV max in suboptimal 
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Fig. 2. Average system spectral efficiency (bit/s/Hz) versus maximum transmit 

'I 

power allowance, P ma x. for different levels of INR, — <f-. 



algorithm 1 is set to 5. Note that if the transmitter is unable 
to guarantee the minimum required power transfer P^', we 
set the system capacity for that channel realization to zero to 
account for the corresponding failure. For the sake of illustra- 
tion, we define the interference-to-signal processing noise ratio 

(INR) as 

A. Average System Spectral Efficiency 

Figure |2] plots the average system spectral efficiency 
(bit/s/Hz) of the three proposed algorithms versus the maximum 

transmit power allowance, P max , for different INRs, It 
can be observed that the spectral efficiency of all algorithms 
increases with increasing P max . This is attributed to the fact that 
the transmitters in the algorithms radiate all the power available 
at the transmitter at every time instant whenever it is possible. 
Besides, an increasing amount of INR impairs the spectral 
efficiency of the system, despite the fact that part of the energy 
can be harvested by the receiver for satisfying the minimum re- 
quired power transfer. On the other hand, suboptimal algorithm 
1 performs very close to the optimal algorithm in all considered 
scenarios due to the iterative optimization. Besides, as expected, 
suboptimal algorithm 2 gives an excellent performance when 
-Pmax is large since it is asymptotically optimal in the high 
SINR regime. Yet, the performance of suboptimal algorithm 2 
is less appealing in the low SINR regime, e.g. for INR = 20 dB 
and P max < 18 dBm, compared to the other two algorithms, 
due to the approximation of the objective function, cf. (fl4l . 

B. Average Harvested Power and Power Splitting Ratio 

Figures [3] and |4] show, respectively, the average power split- 
ting ratio, pi, of the optimal algorithm and the two suboptimal 
algorithms versus maximum allowed transmit power, P max , for 

different levels of INR, For the optimal algorithm, it can be 
observed in Figure [3]that in all considered scenarios, the values 
of pi increase w.r.t. to the maximum transmit power allowance 
Pmax- Although splitting more power for information decoding 
could also possibly increase the associated interference power, 
Pi{& 2 z a + cf ), the power gain due to an increasing pi in the 
signal strength of the desired signal, piPi\Hi\ 2 , is able to 
counteract the performance impairment since the desired signal 
strength dominates the total received power. Besides, the slope 
of the curves depends heavily on the values of INRs and P ma x- 




Fig. 3. Average optimal power splitting ratio, pi, of the optimal algorithm 

versus the maximum transmit power allowance, P ma x, for different levels of 
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Indeed, the role of p t is to balance the spectral efficiency caused 
by piPi \Hi | 2 and the degradation caused by pi (a^ a +c|. ) which 
results in a non-trivial trade-off between INR, P ma x, and pi, 
cf. (O, (f3]l. For instance, for a fixed small value of P max (e.g. 
Pmax < 18 dBm), the average optimal value of pi is higher for 
INR= 20 dB than for INR = 10 dB. Yet, for a fixed large value 
of Pmax (e.g. Pmax > 18 dBm), the average optimal value of pi 
is higher for INR= 10 dB than for INR = 20 dB. The detailed 
trade-off between these variables will be investigated in future 
work. On the other hand, Figure |4] shows that the values of pi 
for the two suboptimal algorithms increase monotonically w.r.t. 
Pnax- In particular, both suboptimal algorithms have a higher 
preference for splitting more power for information decoding in 
the high INR regime. This is because the suboptimal algorithms 
only split enough power to satisfy the minimum required power 
transfer constraint CI with equality (see Figure |5) and allocate 
the remaining power for information decoding. 

Figure reveals that, for the optimal algorithm, the average 
harvested power is a bell-shape function of the maximum 
transmit power allowance P max . Besides, the amount of har- 
vested power is always larger than the minimum required 
power transfer P m ;„. On the contrary, for the two suboptimal 
algorithms, the receivers only harvests just enough power for 
satisfying the minimum required power Pl^ which limits the 
power gain Pipi in the channel capacity. Indeed, the two 
suboptimal algorithms both achieve a lower average system 
capacity and a lower average harvested power compared to 
the optimal algorithm. This is because the two suboptimal 
algorithms are not able to fully exploit the spectral efficiency 
gain achieved by piPi\Hi\ 2 and to reduce the degradation 
caused by Pi{<T 2 a + erf ). 

VI. Conclusions 

In this paper, we formulated the power allocation algorithm 
design for simultaneous wireless information and power transfer 
in OFDM systems as a nonconvex optimization problem. The 
problem formulation took into account the minimum required 
power transfer for system capacity maximization. The prob- 
lem was solved by a one-dimensional full search and convex 
optimization techniques which incurred a high complexity at 
the transmitter. Hence, two low-complexity suboptimal iterative 
algorithms were proposed to find a good compromise between 
computational complexity and system performance. Simula- 




Fig. 5. Average harvested power (dBm) versus the maximum transmit power 
allowance, P max , for different levels of INR, 

a s 

tion results illustrated that the proposed suboptimal algorithms 
achieved a close-to-optimal system capacity. 
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